Electron Spins in Artificial Atoms and Molecules for Quantum Computing * 
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Achieving control over the electron spin in quantum dots (artificial atoms) or real atoms promises 
access to new technologies in conventional and in quantum information processing. Here we review 
our proposal for quantum computing with spins of electrons confined to quantum dots. We discuss 
the basic requirements for implementing spin-qubits, and describe a complete set of quantum gates 
for single- and two-qubit operations. We show how a quantum dot attached to leads can be used for 
spin filtering and spin read-out, and as a spin-memory device. Finally, we focus on the experimental 
characterization of the quantum dot systems, and discuss transport properties of a double-dot and 
show how Kondo correlations can be used to measure the Heisenberg exchange interaction between 
the spins of two dots. 



I. INTRODUCTION 

Coulomb blockade phenomena have atracted much in- 
terest during the last few dacades 0. Creation of con- 
fined electron systems at the nanometer scale has made it 
possible to study the quantum-meachnical nature of the 
band electron in a variety of materials. The tunability 
of the quantum-dot devices provided a unique opportu- 
nity to study the charging effects, and hence, the correla- 
tion effects associated with the Coulomb charging energy. 
However, in confined systems of smaller sizes, where the 
size-quantization energy is resolved, new correlations set 
in due to the Pauli exclusion principle. This brings along 
the electron spin as a degree of freedom in quantum con- 
fined structures, and accessing it in a deterministic way 
would allow for novel implementations in quantum infor- 
mation processing. An increasing number of spin-related 
experiments indeed show that the spin of the elec- 
tron in quantum-confined nanostructures is a promising 
candidate for information processing, due to the unusu- 
ally long (100's nanosecs) spin dephasing times On 
the other hand there are propsed methods || for efficient 
and deterministic control of the spin state in single quan- 
tum dots as well as methods of entangling the spins of 
two dots, the latter being a crucial element in quantum 
information processing. Thus, the field of interest falls 
into two parts, one being improving the technologies for 
conventional computation, and the other - implement- 
ing fundamentally new algorithms of computation with 
quantum bits of information (qubits) and devising a scal- 
able quantum computer in the long run. In conventional 
computers, the electron spin can be expected to enhance 
the performance of quantum electronic devices, such as 
spin-transistors (based on spin-currents and spin injec- 
tion), non- volatile memories, single spin as the ultimate 
limit of information storage etc. (^^]. For implementing 



quantum computing pQ] , as first pointed out in Ref. [g], 
the spin of a confined electron appears as the most nat- 
ural candidate for the qubit. Indeed, provided the spin- 
orbit coupling is negligible, the intrinsic two-state space 
of the spin encodes exactly one qubit and allows for no 
undesired parts of the Hilbert space, transitions to which 
could lead to leakage errors in quantum computation. 
We have shown M that the spin qubits, when located 
in quantum-confined structures such as semiconductor 
quantum dots or atoms or molecules, satisfy all require- 
ments needed for a scalable quantum computer. 

The long distances, of up to 100 /.im ||, over which 
spins can be transported phase-coherently, make the elec- 
tron spin a plausible candidate for quantum informa- 
tion transmission in solid state devices. A spin-qubit 
attached to a mobile electron can be transported along 
conducting wires between different subunits in a quan- 
tum network 11 l^j. Entangled electrons, which can 
be created in coupled quantum dots or via a supercon- 
ductor |D| , provide a source of Einstein-Podolsky-Rosen 
(EPR) pairs Jll];[f^] , wich are necessary for secure quan- 
tum communication. 

The methods used to implement the electron spin in 
conventional computers and in quantum computers are 
often identical, because of the quantum-mechanical na- 
ture of the electron and its spin. Our short-term goal is 
to find ways to control the coherent dynamics of electron 
spins in quantum-confined nanostructures. The use of 
solid state physics as a base for implementing the quan- 
tum computer is motivated by the unparalleled flexibility 
in designing an appropriate medium for the realization of 
a given physical phenomena. 

In the following, we review the status of our theoretical 
efforts towards the goal of implementing quantum com- 
putation with electron spins in quantum-confined nanos- 
tructures. 



Tnvited review prepared for Special Issue of Semiconductor Science and Technology, "Semiconductor Spintronics" , ed. H. 
Ohno, 2002. 
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A. Quantum Computing and Quantum Dots 

The possibility of outperforming classical computation, 
which opens up in quantum algorithms such as the one 
discovered by Shor H] and by Grover |15|, has attracted 
much interest. A quantum algorithm makes use of the 
quantum computers's ability to exist in any superposition 
of the states of its binary basis and to perform quantum 
time evolution for computation; hence the parallelism of 
quantum computing. The requirement for the quantum 
bit of information (qubit), which is at the heart of the 
quantum computer, is that it can exist in any state of 
a quantum two two-level system, i.e. \ip) — a\0) + P\l), 
where |0) and |1) are the states of the "classical" bit, and 
|cv| 2 + |/3| 2 = 1. Apart from this, a qubit should be able to 
couple to any other qubit in the quantum comupter and 
form a coherent two-qubit state. These two elements, 



which are the one- and two-qubit gates, are sufficient for 
forming a many-qubit coherent state and implementing 
any quantum algorithm. 

A recently growing list of quantum tasks [|ll],[l6| such 
as cryptography, error correcting schemes, quantum tclc- 
portation, etc. have indicated even more the desirability 
of experimental implementations of quantum computing. 
On the other hand, there is also a growing number of 
proposed physical implementations of qubits and quan- 
tum gates. A few examples are: Trapped ions JlT 
ity QED " 
vices lEll 



nuclear spins |19|,|20[, superconducting de- 
and our qubit proposal Q based on the 
spin of the electron in quantum-confined nanostructures, 
and in particular in quantum dots with an all-electrical 
control of spin. Subsequent proposals such as pp|,p5| , p6| 




are based on the same principles as introduced in 
reviewed herein. 



and 




high-g layer quantum well 

FIG. 1. Quantum dot array, controlled by electrical gating. The electrodes (dark gray) define quantum dots (circles) by 
confining electrons. The spin 1/2 ground state (arrow) of the dot represents the qubit. These electrons can be moved by 
electrical gating into the magnetized or high-g layer, producing locally different Zeeman splittings. Alternatively, magnetic 
field gradients can be applied, as e.g. produced by a current wire (indicated on the left of the dot-array). Then, since every 
dot-spin is subjected to a different Zeeman splitting, the spins can be addressed individually, e.g. through ESR pulses of an 
additional in-plane magnetic ac field with the corresponding Larmor frequency ojl = g/J.BB±/h. Such mechanisms can be used 
for single-spin rotations and the initialization step. The exchange coupling between the quantum dots can be controlled by 
lowering the tunnel barrier between the dots. In this figure, the two rightmost dots are drawn schematically as tunnel-coupled. 
Such an exchange mechanism can be used for the XOR gate operation involving two nearest neighbor qubits. The XOR 
operation between distant qubits is achieved by swapping (via exchange) the qubits first to a nearest neighbor position. The 
read-out of the spin state can be achieved via spin-dependent tunneling and SET devices , or via a transport current passing 
the dot J27]]. Note that all spin operations, single and two spin operations, and spin read-out, are controlled electrically via the 
charge of the electron and not via the magnetic moment of the spin. Thus, no control of local magnetic fields is required, and 
the spin is only used for storing the information. This spin-to-charge conversion is based on the Pauli principle and Coulomb 
interaction and allows for very fast switching times (typically picoseconds) . A further advantage of this all-electrical scheme is 
its scalability into an array of arbitrary size. 
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Semiconductor quantum dots are structures where 
charge carriers are confined in all three spatial dimen- 
sions, the dot size being of the order of the Fermi wave- 
length in the host material, typically between 10 nm and 
1 /im Q. The confinement is usually achieved by elec- 
trical gating of a two-dimensional electron gas (2DEG), 
possibly combined with etching techniques, see Fig. [fl. 
Precise control of the number of electrons in the conduc- 
tion band of a quantum dot (starting from zero) has been 
achieved in GaAs heterostructures (28). The electronic 
spectrum of typical quantum dots can vary strongly when 
an external magnetic field is applied |l],^8| , since the mag- 
netic length corresponding to typical laboratory fields 
B w 1 T is comparable to typical dot sizes. In coupled 
quantum dots Coulomb blockade effects J2{|, tunneling 
between neighboring dots p|p9[, and magnetization |3(J] 
have been observed as well as the formation of a delocal- 
ized single-particle state |}l]] . 



II. GENERAL CONSIDERATIONS FOR 
QUANTUM COMPUTING WITH SPINS 

A. Coherence 

Magneto-optical experiments, based on time-resolved 
Faraday rotation measurements, show long spin coher- 
ence times in doped GaAs in the bulk and a 2DEG ||. 
At B = and T = 5 K, a transverse spin lifetime (de- 
phasing time) T 2 * exceeding 100 ns was measured, with 
experimental indications that this time is a single-spin ef- 
fect ||. Since this number still includes inhomogeneous 
effects — e.g. g-factor variations in the material, leading 
to spins rotating with slightly different frequencies and 
thus reducing the total magnetization — it represents only 
a lower bound of the decoherence time T2 of a single 
spin, Ti > T 2 *, which is relevant for using spins as qubits. 
Using the same pump-probe technique, spin dephasing 
times in semiconductor (CdSe) quantum dots have been 
measured [p2| , with at most one spin per dot. The rel- 
atively small T 2 * dephasing times (a few ns at vanishing 
magnetic field), which have been seen in these experi- 
ments, probably originate from a large inhomogeneous 
broadening due to a strong variation of g-factors fl32"| . 
Nevertheless, the fact that many coherent oscillations 
were observed |32 provides strong experimental support 
to the idea of using electron spin as a qubit. 



B. Upscaling 

To outperform a classical computer, a quantum com- 
puter will need a number of qubits on the order of 10 5 . 
Hence, it is essential that the underlying concept can 
be scaled up to a large number of qubits. This scal- 
ing requirement is, in principle, achievable with spin- 
based qubits confined in quantum dots, since producing 



arrays of quantum dots |ll|]33j ] is feasible with today's 
technology of defining nanostructures in semiconductors. 
Of course, the actual implementation of such arrays (see 
Fig. ^) including all the needed circuits poses tremendous 
experimental challenges, but at least we are not aware of 
any physical restriction which would exclude such an up- 
scaling for spin-qubits. 



C. Switching 

Quantum gate operations ca n be controlled through an 
effective Hamiltonian (see Sec. pill and IV ) 



H{t) = M*) S * ' S o +^^B9i{t) B<(t) ■ 8i . (1) 

i<j i 

The coupling constants Jy , gi and the magnetic field (lo- 
cal) Bi are controlled via external gate fields, which are 
switched with some pulses v(t). In the following we as- 
sume Jij to be non-zero only for the neighboring qubits. 
Note, however, that in cavity-QED systems there is also 
a long-range coupling of qubits H] , and that long-range 
coupling via a superconductor is also possible [^5|. But 
even if the exchange coupling is only local, operations on 
non-neighboring qubits can still be performed. This is 
achieved by swapping states of neighboring qubits (see 
Sec. IV), which allows one to move the qubit around in 
an array of quantum dots and couple it to the desired 
other qubit. 



For the gating mechanisms described in Sec. Ill and IV, 
only the time integral J Q P(v(t))dt (mod 2ir) is impor- 
tant. Here, P(v(t)) stands for the exchange coupling J 
or the Zeeman interaction. The requirement on the pulse 
shape is that it does not violate the validity of the effec- 
tive Hamiltonian (|l|), but otherwise the gating mecha- 
nisms are independent of the actual shape of v(t). Since 
the effective Hamiltonian (|l|) was obtained by project- 
ing out h igher energy states of one and two coupled dots 
(see Sec. IV A), care should be taken that the pulses do 
not excite the quantum dots to the projected-out energy 
levels. This can be achieved by switching v(t) adiabat- 
ically, i.e. such that \v/v\ <C Se/h, where Se is the en- 
ergy scale on which excitations may occur. We find that 
v(t) = vq sech(t/Ai), where vq is the pulse amplitude 
and At the characteristic width, is optimal for a fast adi- 
abatic switching, provided 1/At <C 5elh. For a detailed 
analysis of adiabatic switching, see [|36| . 

A single qubit operations can be performed for ex- 
ample in g-factor-modulated materials, as described in 
Sec. [II. A spin can be rotated by a relative angle of 
(p. = Ag c HfiBBT /2h through changing the effective g- 
factor by Ag c g for a time r. Thus, a typical switching 
time for an angle ip = 7r/2, a field B = 1 T, and Ag e g ss 1 
is t s w 30 ps. If slower operations are required, they are 
easily implemented by choosing a smaller Ag ff, reduc- 
ing the magnitude of the field B, or by replacing tp by 
ip + 2im with integer n. 
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Next we consider two exchange-coupled spins, which 
perform a square-root-of-swap gate for the integrated 
pulse J Q S J(t)dt/h = 7r/2, as described in Sec. [V. We 
apply a pulse J{t) = Josech((i — r s /2)/Ai) with Jo = 
80 /ieV, and choose At — 4 ps, which gives for the 
switching time t s ps 30ps, and the adiabaticity criterion 
h/Atn 150^eV< Se. 



D. Error Correction 

Realization of a reliable error-correction scheme ]37| ] 
is one of the main goals in quantum computation. The 
known schemes for fault-tolerant quantum computation 
work if the gate operation error rate does not exceed a 
certain threshold value, usually about 10 -4 (depending 
on the scheme) [ 38|. If we take the ratio of the switching 
times fro m Sec . [I Q , t s ss 30 ps, and the dephasing time 
from Sec. 



II A 



100 ns, we obtain a value close to 
this threshold. Thus, in our proposal, we can expect an 
arbitrary upscaling of the quantum computer, and we are 
no further limited by decoherence and lacking gate pre- 
cision. We note that implementing an error-correction 
scheme requires a larger number of gate operations, and 
therefore, it is desirable to perform them in parallel; oth- 
erwise the pursued gain in computational power is used 
up for error correction. Hence, one favors concepts where 
a localized control of the gates can be realized such that 
operations can be performed in parallel. However, since 
there are still many milestones to reach before sophisti- 
cated error-correction schemes can be applied, one should 
by no means disregard setups where gate operations are 
performed in a serial way. 



III. SINGLE-SPIN ROTATIONS 



For quantum computing it is necessary (but not suffi- 
cient) to perform one-qubit operations. In the context of 
spin-qubits, it translates into single-spin rotations. This 
can be achieved by exposing a specific qubit to a time- 
varying Zeeman coupling (g/isS-B)(t) ]38[ j , which can be 
controlled through both the magnetic field B and/or the 
g-factor g. Since only phases have a relevance, it is suffi- 
cient to rotate all spins of the system at once (e.g. by an 
external field B), but with a different Larmor frequency. 

Localized magnetic fields can be generated with the 
magnetic tip of a scanning force microscope, a magnetic 
disk writing head, by placing the dots above a grid of 
current-carrying wires, or by placing a small wire coil 
above the dot etc. 

Single-spin rotations can be achieved by ESR tech- 
niques p9| ]. One applies a static local magnetic field B 
for the qubit(s), which should be rotated. An ac mag- 
netic field is then applied perpendicular to the first field 
with the resonant frequency that matches the Larmor 



frequency lul — g^iBB/h. Due to paramagnetic reso- 
nance pOjl , this causes spin- flips in the quantum dots with 
the corresponding Zeeman splitting. 

The equilibrium position of the electron can be moved 
around through electrical gating. Thus, if the electron 
wave function is pushed into a region with a different 
magnetic field strength or (effective) g-factor, one pro- 
duces a relative rotation around the direction of B by an 
angle of tp = (g'B' — gB)(iBi~/2h, see Fig. |l]. Regions 
with an increased magnetic field can be provided by a 
magnetic (dot) material while an effective magnetic field 
can be produced e.g. with dynamically polarized nuclear 
spins (Overhauser effect) [ |39[ . 

We shall now explain a concept for using g-factor- 
modulated materials jll],[33| . In bulk semiconductors the 
free-electron value of the Lande g-factor go = 2.0023 is 
modified by spin-orbit coupling. Similarly, the g-factor 
can be drastically enhanced by doping the semiconduc- 
tor with magnetic impurities |)|,[lj . In confined structures 
such as quantum wells, wires, and dots, the g-factor is 
further modified and becomes sensitive to an external 
bias voltage pi] ]. We have numerically analyzed a sys- 
tem with a layered structure (AlGaAs-GaAs-InAlGaAs- 
AlGaAs), in which the effective g-factor of electrons is 
varied by shifting their equilibrium position from one 
layer to another by electrical gating. We have found that 
in this structure the effective g-factor can be changed by 
about Ag e s « 1 [Q. Such a gate-controlled g-factor 
modulation has now been confirmed experimentally p2| . 



IV. TWO-QUBIT GATES 

The main component for every computer concept is a 
multi-(qu)bit gate, which eventually allows calculations 
through combination of several (qu)bits. Since two-qubit 
gates are (in combination with single-qubit operations) 
sufficient for quantum computation p3[ — they form a 
universal set — we now focus on a mechanism that cou- 
ples pairs of spin-qubits. Such a mechanism exists in 
coupled quantum dots, resulting from the combined ac- 
tion of the Coulomb interaction and the Pauli exclusion 
principle. Two coupled electrons in absence of a magnetic 
field have a spin-singlet ground state, while the first ex- 
cited state in the presence of strong Coulomb repulsion is 
a spin triplet. Higher excited states are separated from 
these two lowest states by an energy gap, given either 
by the Coulomb repulsion or the single-particle confine- 
ment. The low-energy dynamics of such a system can be 
described by the effective Heisenberg spin Hamiltonian 



H B {t) = J{t) Si ■ S 2 , 



(2) 



where J(i) denotes the exchange coupling between the 
two spins Si and S2, i-e. the energy difference between 
the triplet and the singlet. After a pulse of J it) with 
/J*" dtJ(t)/h = Jot s /Tl = 7r (mod 2n), the time evolu- 
tion U(t) = Texp(i J Q H s (T)d,T/h) corresponds to the 
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"swap" operator U sw , whose application leads to an in- 
terchange of the states in qubit 1 and 2 || . While U sw is 
not sufficient for quantum computation, any of its square 

roots Usl 2 , say uU 2 \4>x) = (\4>x) + i\x<t>)) I '(1 + i), turns 
out to be a universal quantum gate. Thus, it can be 
used, together with single-qubit rotations, to assemble 
any quantum algorithm. This is shown by constructing 
the known universal gate XOR , through combination 

of uU, 2 and single-qubit operations exp(iTtSf /2), applied 
in the sequence H, 

Uxor = ^ 12)31 e-^ s i U^ 2 U^ 2 - (3) 



With these universal gates at hand, we can reduce 
the study of general quantum computation to the study 
of single-spin rotations (see Sec. DTD) and the exchange 
mechanism, in particular how J(t) can be controlled 
experimentally. The central idea is that J(t) can be 
switched by raising or lowering the tunneling barrier be- 
tween the dots. In the following, we shall review our 
detailed calculations to describe such a mechanism. We 
note that the same principles can also be applied to 
other spin systems in quantum-confined structures, such 
as coupled atoms in a crystal, supramolecular structures, 
and overlapping shallow donors in semiconductors [ p0|j26| 
etc., using similar methods as explained below. 



A. Coupled Quantum Dots 

We consider a system of two tunnel-coupled quantum 
dots, achieved by gating a two-dimensional electron gas 
(2DEG), as described in Sec. [A. The dots are arranged 
in a plane (see Fig. ||), at a sufficiently small distance 2a, 
such that the electrons can tunnel between the dots. The 
tunnel junction between the dots, as well as the number 
of electrons on each dot, is controlled by the depicting 
gates (see Fig. |l|). We consider the case of similar dots, 
and each dot contains an odd number of electrons with 
a spin 1/2 ground state. Furthermore, we simplify our 
consideration by retaining only one electron per dot and 
assuming that the rest of the electrons form a closed shell 
and merely contribute to the confining potential of the 
gates. We model the two coupled dots with the Hamil- 

where 



tonian |39| H = £ i=1 _ 2 hi + C + H z = H olh + H z , 
the single-electron dynamics in the 2DEG (xy-pl&ne) is 
described through 



1 / e 
h = — Pi - -A(rj) 

Am \ c 



(4) 



with m being the effective mass and W(jCi) the con- 
finement potential as given below. A magnetic field 
B = (0, 0, B) is applied along the z-axis, which cou- 
ples to the electron spin through the Zeeman interac- 
tion Hz and to the charge through the vector potential 
A(r) = (— y, x, 0). In almost depleted regions, like few- 
electron quantum dots, the screening length A can be ex- 
pected to be much larger than the screening length in 



bulk 2DEG regions (where it is 40 nm for GaAs). Thus, 
for small quantum dots, say A 3> 2a « 40 nm, we need to 
consider the bare Coulomb interaction C — e 2 / «|ri — r2 1, 
where n is the static dielectric constant. The confinement 
and tunnel- coupling in Eq. (|J) for laterally aligned dots 
is modeled by the quartic potential 



W(x,y) = 



1 
4a~ 2 



a 2 ) 2 + y 2 



(5) 



with the inter-dot distance 2a and as = yJh/muiQ the ef- 
fective Bohr radius of the dot. Separated dots (a > as) 
are thus modeled as two harmonic wells with frequency 
lu . This is motivated by the experimental evidence that 
the low-energy spectrum of single dots is well described 
by a parabolic confinement potential Eq]. 




FIG. 2. Double dot geometry 

Now we first classify the two-particle states accord- 
ing to the available symmetries, and then we calculate 
the phenomenological parameters within our toy model. 
The wave function of the two electrons in the double dot 
(DD) can be chosen to be a product of an orbital and 
a spin part. Then, in a singlet (triplet) state the or- 
bital part is symmetric (antisymmetric) with respect to 
the interchange of the electrons, while the spin part is 
antisymmetric (symmetric). With respect to the mir- 
ror reflection in the yz-plane [ ffo| (see Fig. ^) the single- 
particle states fall into two symmetries, which we label 
by the quantum number n = ±. The energy difference 
between the state with n = — and that with n = + is 
given by 2to > (to being the interdot hopping ampli- 
tude). Within the low energy sector, the lowest singlet 
state and the triplet states are then given by [M 



|00) = 



1 



|ii) = dL T < T |o) , |i-i) = dL i 4ilo> 

|10) = -^(dL T <;+dL;< T )|0) , 



(6) 



where the notation \SS Z ) stands for the angular momen- 
tum representation of the total spin of the two electrons. 
Here, the second quantized operator creates an elec- 
tron in the orbital state n with spin a. The vacuum 
state |0) includes the disregarded electrons. The interac- 
tion parameter (f> depends on the interplay between the 
tunneling and Coulomb interaction. We have calculated 
6 Wh within the Hund-Mulliken method, 
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(7) 



where i# and J/h are the extended inter-dot tunnel- 
ing amplitude and on-site Coulomb repulsion, respec- 
tively 1 39 1 . We note that tji = to + tc, with the contri- 
bution tc coming from the Coulomb interaction J39j ] and 
vanishing with vanishing to. For detached dots we have 
cj) = 1, and tfi < 1 occurs due to double occupancies in 
the dots, and 1 —> for vanishing Coulomb interaction. 
By varying a, we plot cj) versus to on Fig. ||. Next, we 
note that the singlet state in (^) represents an entangled 
state of two electrons, with the entanglement being 



V = 



20 



1 



(8) 



Here, we used the measure of entanglement introduced 
in Ref. M. 
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FIG. 3. Parameter cj> versus the interdot hopping ampli- 
tude to at different magnetic fields (Hund-Mulliken calcula- 
tion). We used GaAs quantum dots with confinement energy 
Tiojq = 3meV and dielectric constant k = 13.1. 

Next, assuming adiabatic switching of the coupling 
constants, we arrive at the effective Hamiltonian (g) by 
means of the following mapping p8| , fl7[ 

^2 dla cr crcr>dn><j> = S+Snn' + ( "7T S - +«n0_Tj 8- nn >, 



Sr Sa- 



n Sj ± S 2 , T = Sj x S 2 , <j> ± = 
b 2 ), and <r are the Pauli matrices. 



= V2(l ± 
Mapping 



where S± 

@ projects out the higher energy sector of the DD and 
keeps only the states (|J). We note that the spin 1/2 
operators Si^ are nothing but the intermixed electron 



spins, and hence represent the spin degrees of freedom 
of the DD. When detaching the two dots adiabatically, 
one always obtains one electron in each dot |56| and the 
spins Si^ then stand for the true electron spins. Also 
note that during an adiabatic coupling of the two qubits 
(dots) each of the spins Si^ carries the initial informa- 
tion of its qubit, which gets modified only through the 
Heisenberg exchange interaction J , see Eq. (Q) . The adia- 
baticity criterion discussed in Sec. II C applies hence here 
with Se — min(?kJo, Uh). 

The Heisenberg exchange interaction J, which is de- 
fined as the energy difference between the triplet and 
singlet states (^), is the only parameter of interest for 
the two-qubit dynamics (provided the adiabaticity cri- 
terion is fulfilled). Although for a real structure it is 
best to have methods to measure J for different values 
of the gate voltages, i.e. to characterize the structure 



experimentally (see Sec. VI), we still find it instructive 
to analyse various contributions to J within our realistic 
model of the DD. In particular, we show that breaking 
the time-reversal symmetry by means of applying a mag- 
netic field leads to a singlet-triplet transition in the DD. 
We calculate J using different methods and compare the 
results. A generic expression for J is straightforwardly 
obtained from the states (|6|), 



J= V 



1 - 



1 



) 2 U 



H 



1 



(10) 



with V = (T\C\T) - (00+|C|00+), t H = t - 
(00+|C|0(T)/2, and U H = (0tr|C|00-) - (00+|C|00+). 
Here, |T) stands for any of the triplet states in (|J), and 
|00 ± ) denotes the singlet state of (^ taken at (b = ±1. 
In the Hund-Mulliken approach, <p is given by (0) and 
expressions for V, tn, and Uh were obtained in Ref. [ p9[ . 
The Heisenberg exchange interaction then reads 



J = V 



2 



(ii) 



We note that the component V is responsible for mak- 
ing J ferromagnetic, i.e. J < 0. In the standard Hub- 
bard approach for short-range Coulomb interaction it is 
assumed that there is no overlap between the electron 
wave function on different dots, though there is a finite 
hopping amplitude t (chain model) . The exchange terms 
of the Coulomb interaction hence vanish, and this cor- 
responds to setting V — > 0, tn — > t, and Uh — > U, 
with U being the on-site (short-range) Coulomb repul- 
sion | j39[ . The parameter <p is given by a formula analo- 
gous to (Q) . The Heisenberg exchange interaction is then 
always antiferromagnetic, J = */ {U/2) 2 + 4t 2 — U/2 > 0. 
Finally, in the Heitler-London approach, the double oc- 
cupancy of each dot is neglected, which corresponds to 
setting <fr = (1 — S)/(l + S), where S is the overlap in- 
tegral between the electron wave functions on the two 
dots. Note that the Heitler-London method gives qual- 
itatively wrong results for both the case of strong and 



G 



weak Coulomb interaction, since the interaction param- 
eter (j> is not sensitive to the Coulomb interaction. How- 
ever, for intermediate strengths of the Coulomb interac- 
tion (Ujj ~ hu>o) it gives qualitatively correct results. 
The parameters V, tn, and Uh are identical for both 
Heitler-London and Hund-Mulliken methods. We find it 
appropriate to use the Heitler-London approximation for 
simplicity in a GaAs system. Formula (|l0|) then reduces 

to M, 



J 



hu 



sinh (2d 2 



26-1 N 



(12) 



zVb 



e""" I (bd 2 ) - e d ^ b -^ b I (d 2 



6-1 



where we have introduced the dimensionless distance 
d = a/aB between the dots and the magnetic compression 
factor 6 = B/Bq = \fl + lo 2 l /loq with the Larmor fre- 
quency io l — eB/2mc. The zeroth order Bessel function 
is denoted by I . In Eq. fll^), the first term comes from 
the confinement potential, while the terms proportional 
to the parameter c = \/tt /2(e 2 / kob) /hvo result from the 
Coulomb interaction C; the exchange term is recognized 
by its negative sign. We plot J [Eq. rtjj)] in Fig. f| as a 
function of B and d. We note that J(B = 0) > 0, which 
is generally true for a two-particle system with time- 
reversal invariance. We observe that over a wide range 
of the parameters c and a, the sign of J(B) changes from 
positive to negative at a finite value of B (for the param- 
eters chosen in Fig. ||(a) at B w 1.3 T). J is suppressed 
exponentially either by compression of the electron or- 
bitals through large magnetic fields (6 3> 1), or by large 
distances between the dots (d^> 1), where in both cases 
the orbital overlap of the two dots is reduced. This ex- 
ponential suppression, contained in the 1/sinh prefactor 
in Eq. (|l^), is partly compensated by the exponentially 
growing exchange term oc exp(2d 2 (6 — 1/6)). In total, 
J decays exponentially as exp(— 2<i 2 6) for large b or d. 
Since the sign reversal of J — signalling a singlet-triplet 
crossing — results from the long-range Coulomb interac- 
tion, it is not contained in the standard Hubbard model 
which takes only short-range interaction into account. 
Fig. |4| compares the results of different methods. 

We remark again that the exponential suppression of 
J is very desirable for minimizing gate errors. In the ab- 
sence of tunneling between the dots we still might have 
direct Coulomb interaction left between the electrons. 
However, this has no effect on the spins (qubit) pro- 
vided the spin-orbit coupling is sufficiently small, which 
is the case for s-wave electrons in GaAs structures with 
unbroken inversion symmetry (this would not be so for 
hole-doped systems since the hole has a much stronger 
spin-orbit coupling due to its p-wave character). For a 
detailed discussion of spin-orbit interaction and its effect 
on exchange and XOR gates we refer to Ref . M . 
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FIG. 4. Exchange coupling J (full line) for GaAs quan- 
tum dots with confinement energy Tioj = 3 meV and c = 2.42. 
For comparison we plot the usual short-range Hubbard result 
J = At 2 /U (dashed-dotted line) and the extended Hubbard 
result |39| J — V + 4t 2 H /UH (dashed line). In (a), J is plot- 
ted as a function of the magnetic field B at fixed inter-dot 
distance d — a/as = 0.7, while in (b) as a function of the 
inter-dot distance d = a/as at B = 0. 



V. MEASURING A SINGLE SPIN (READ-OUT) 

A. Spin Measurements through Spontaneous 
Magnetization 

One scheme for reading out the spin of an electron 
on a quantum dot is implemented by tunneling of this 
electron into a supercooled paramagnetic dot M . There 
the spin induces a magnetization nuclcation from the 
paramagnetic metastablc phase into a ferromagnetic do- 
main, whose magnetization direction (8, tp) is along the 
measured spin direction and which can be measured by 
conventional means. Since this direction is continuous 
rather than only one of two values, we describe this 
generalized measurement in the formalism of positive- 
operator-valued (POV) measurements |5(]] as projection 
into the overcomplete set of spin- 1/2 coherent states 
\6,ip) = cos(0/2)||) + e^sin(6'/2)||). Thus if we in- 
terpret a magnetization direction in the upper hemi- 
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sphere as | f), we have a 75%- reliable measurement, since 
(1/2tt) J g>7r ^ 2 dQ\(^\9 7 ip)\ 2 = 3/4, using the normaliza- 
tion constant 2ir for the coherent spin states. 



B. Quantum Dot as Spin Filter and 
Read-Out/Memory Device 

We discuss now a setup — quantum dot attached to in- 
and outgoing current leads I = 1, 2 — which can be oper- 
ated as a spin filter, or as a read-out device, or as a spin- 
memory where a single spin stores the information [H . 

A new feature of this proposal is that the spin- 
degeneracy is lifted with different Zeeman splittings in 
the dot and in the leads, e.g. by using materials with dif- 
ferent effective g- factors for leads and dot ||] . This results 
in Coulomb blockade peaks and spin-polarized currents 
which are uniquely associated with the spin state on the 
dot. 

The setup is described by a standard tunneling Hamil- 
tonian H a + Ht [|5l]], where H a = Hl + Hr> describes 
the leads and the dot. Ho includes the charging and in- 
teraction energies of the electrons in the dot as well as 
their Zeeman energy ±gfi B B/2 in an external magnetic 
field B. The tunneling between leads and the dot is de- 
scribed by H T = J2i,k, P , s tipCiksdps + h.c, where ci ks an- 
nihilates electrons with spin s and momentum k in lead I, 
and d ps annihilates electrons in the dot. We consider the 
Coulomb blockade regime jij where the charge on the dot 
is quantized. Then we apply a standard master-equation 
approach |52p|] with a reduced density matrix of the dot 
and calculate the transition rates in a "golden-rule" ap- 
proach up to 2nd order in Ht- The first-order contri- 
bution to the current is the sequential tunneling current 
I s , where the number of electrons on the dot fluctu- 
ates and thus the processes of an electron tunneling from 
the lead onto the dot and vice versa are allowed by en- 
ergy conservation. The second-order contribution is the 
cotunneling current I c |53j |, involving a virtual interme- 
diate state with a different number of electrons on the 
dot. 

We now consider a system, where the Zeeman splitting 
in the leads is negligible (i.e. much smaller than the Fermi 
energy), while on the dot it is given as A 2 = [i B \gB\. We 
assume a small bias A/i = /ii — /i2>0 between the leads 
at chemical potential /ii, 2 and low temperatures so that 
A/i, kT < 8, where S is the characteristic energy-level 
distance on the dot. First we consider a quantum dot in 
the ground state, filled with an odd number of electrons 
with total spin 1/2, which we assume to be 1 1 ) and to 
have energy — 0. If an electron tunnels from the lead 
onto the dot, a spin singlet is formed with energy Es, 
while the spin triplets are (usually) excited states with 
energies Et ± and Et - At the sequential tunneling res- 
onance, fix > Es > /i2, where the number of electrons 
on the dot fluctuates between N and N + 1, and in the 
regime Et — Es, A z > A/i, kT, energy conservation 



only allows ground state transitions. Thus, spin-up elec- 
trons are not allowed to tunnel from lead 1 via the dot 
into lead 2, since this would involve virtual states |T + ) 
and ||), and so we have / s (t) = for sequential tun- 
neling. However, spin down electrons may pass through 
the dot in the process l(T)i — * (Q)/, followed by (fj); — > 
(T) j/. Here the state of the quantum dot is drawn inside 
the circle, while the states in the leads are drawn to the 
left and right, resp., of the circle. This leads to a spin- 
polarized sequential tunneling current / s = I S {1), which 
we have calculated as || 

I s {l)/h=0{ l ix-Es)-e{n 2 -E s ), k B T<A^, (13) 



IeiD/Io 



^cosh- 2 



2k B T 



k B T > A/i, (14) 



where fi = (/i X + /t2)/2 and I = e^i^ 2 /(7i + 12}- Here 
7/ = 2i^v\Ai nn i | is the tunneling rate between lead I and 
the dot, and we have introduced the matrix elements 
Ain'n = Ytpstipin'ldpA 71 )- Similarly, for N even we find 
Is (I) = 0, while for / s (f) a similar result holds |9| as in 
Eqs. @, @. 

Even though I s is completely spin-polarized, a leakage 
of current with opposite polarization arises through co- 
tunneling processes ; still the leakage is small, and the 
efficiency for A z < \Et + — Es\ for spin filtering in the 
sequential regime becomes H (for 71 ~ 72) 



(71 + 72) max{fcsT, A/i} ' 



(15) 



and equivalently for I s (1) / I c (l) at the even-to-odd tran- 
sition. In the sequential regime we have 7^ < k B T,Afi, 
thus, for k B T, A/i < A z , we see that the spin-filtering is 
very efficient. 

We discuss now the opposite case where the leads are 
fully spin polarized with a much smaller Zeeman split- 
ting on the dot [|| . Such a situation can be realized with 
magnetic semiconductors (with effective g-factors reach- 
ing 100 Q) where spin- injection into GaAs has recently 
been demonstrated for the first time . Another possi- 
bility would be to work in the quantum Hall regime where 
spin-polarized edge states are coupled to a quantum dot 
p4j . In this setup the device can be used as read-out for 
the spin state on the dot. Assume now that the spin po- 
larization in both leads is up, and the ground state of the 
dot contains an odd number of electrons with total spin 
1/2. Now the leads can provide and absorb only spin-up 
electrons. Thus, a sequential tunneling current will only 
be possible if the dot state is | J.) (to form a singlet with 
the incoming electron, whereas the triplet is excluded by 
energy conservation). Hence, the current is much larger 
for the spin on the dot being in ||) than it is for ||). 
Again, there is a small cotunneling leakage current for 
the dot-state 1 1), with a ratio of the two currents given 
by Eq. ( plf ) with A z replaced by Et+ — Es- Thus, we 
can probe (read out) the spin-state on the quantum dot 
by measuring the current which passes through the dot. 
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Given that the sequential tunneling current is typically 
on the order of 0.1 — 1 nA |Q, we can estimate the read- 
out frequency I /2ire to be on the order of 0.1 — 1 GHz. 
Combining this with the initialization and read-in tech- 
niques, i.e. ESR pulses to switch the spin state, we have 
a spin memory at the ultimate single-spin limit, whose 
relaxation time is just the spin relaxation time T\. This 
relaxation time can be expected to be on the order of 
100's of nanoseconds jj), and can be directly measured 
via the currents when they switch from high to low due to 
a spin flip on the dot Q . Furthermore, the spin decoher- 
ence time Ti can also be measured via the current, if an 
ESR field is applied to the dot in either sequential tun- 
neling or cotunncling r egim e with normal (unpolarizcd) 
leads, as shown in Ref. 1271. 



VI. ACCESSING THE SINGLET-TRIPLET 
SPLITTING IN DOUBLE DOTS 

Transport measurements can be used to characterize 
a quantum dot system experimentally. A main param- 
eter of interest for quantum computing is the exchange 
interaction J between the spins of two neighboring dots. 
We have considered a setup Q , consisting of two lateral 
quantum dots connected in series between two metallic 
leads, see Fig. || A magnetic field B, applied perpendic- 
ular to the plane of the dots, is used to tune the exchange 
interaction J. A common gate (not shown), with the gate 
voltage V g , can be used to change the electron occupa- 
tion number of the double dot (DD). The conductance 
of the DD versus V g shows peaks of sequential tunneling 
separated by Coulomb blockade valleys [here we consider 
temperatures smaller than the Coulomb correlation en- 
ergy]. We focus on the valley with two electrons in the 
DD. Our consideration also holds for a larger ocupation 
number M = 27V, with N being odd, provided N — 1 
electrons on each dot form a closed shell and can be dis- 
regarded. The two (outer shell) electrons are confined by 
the DD potential, however their spin degrees of freedom 
can be correlated on a much smaller energy scale J. Our 
aim is to provide ways of accessing the exchange inter- 
action J between the spins of the two electrons in the 
DD. 

Readily the differential cotunneling conductance 
through the DD shows distinct features (steps at the bias 
A/i = ± J), which allow one to measure J experimentally. 
However, attaching leads to the DD shifts the energy lev- 
els, and hence, modifies J. Moreover, measuring a small 
value of J requires low temperatures at which the Kondo 
correlations in the leads may be important. We show 
that such Kondo correlations also introduce a correction 
to J, which is temperature dependent. We find that the 
peculiar features in transport properties are better pro- 
nounced in the Kondo regime. For example, the linear 
conductance of the DD as function of temperature shows 
a maximum at temperature T ~ J > 0, which is pro- 
nounced only in the Kondo regime. This maximum can 



be used as an alternative way of measuring J, having 
the advantage that, in the linear regime, the DD is not 
affected by the applied bias. 
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FIG. 5. Double-dot system containing two electrons and 
being coupled in series to two metallic leads at chemical 
potentials pa and pl with bias A/i = pi — pa. The 
electron spins Sl, Sr interact via the exchange interaction 
J = Et — Es, where Et,s is the triplet /singlet energy. 



A. Cotunneling through two tunnel-coupled 
quantum dots 

In the Coulomb blockade regime, the fluctuations of 
the number of electrons on the DD are strongly sup- 
pressed by the Coulomb blockade gap. The conductance 
through the DD is dominated by processes with a vir- 
tual occupation of the DD by a lead electron (hole), i.e. 
cotunncling processes. We consider a realistic DD with 
long range Coulomb interaction between the two elec- 
trons on the DD, as discussed in Sec. [TV A . Using a 
Hund-Mulliken approach, we write down the one-electron 
states as ip±^ = XoiV-a ± f+a)/ ± s ) > where x<j 
is the spinor, ip± a are the lowest orbitals of single dots 
situated at x = ±a, and S — (ifi± a \ l P^a) is the over- 
lap integral. The lowest in energy singlet and triplet are 
given then by Eq. (^]), and the interaction parameter <fi 
is calculated according to (0). The attached leads are 
described by Hi = Y, a ka e kc\ ka c akrTl where c' aka cre- 
ates an electron with momentum k and spin a in lead 
a = L,R. The tunneling between the DD and the leads 



is described by H T = Y naka {t an c akrT d n(J + h.c). Here, 
d ri a annihilates an electron in the state i/w- The tun- 
neling amplitudes are given by tL,± = thj '\/2(l ± S) and 
tR t ± = ±tn/ \/2(l ± S), with t a being the apmlitude to 
tunnel from lead a onto the the adjusting dot at to = 0. 
We map our problem onto a two-level system, with level 1 
corresponding to the singlet state and level 2 to the three 
triplet states. The occupation probabilities of these two 
levels, pi and p2, are given by 



1 



Pi = 



i + 3cxp(-j/r eff ) 



= 1 — P2 



(16) 







The effective temperature T c g depends on the applied 
bias Ap, and thus, describes the heating effects on the 
DD. Solving a master equation for p\ and pi in the co- 
tunneling regime, we find 



1 



1 + A(J/T, Ap/T)^ 2 



= - --In 
T cff T J n 1 + \(-J/T,A(i/T)cj) 2 _ 



(17) 



where A(it, v) = 
the parameter 



sinh(v) tanh(-u/2) — (v/u) tanh(u/2) 



4 l-tanh(u/2) ' 



coshfn') — c oshf v ) 



was introduced in Sec. IV A 



and 
We note 

that the heating effect depends on the interaction param- 
eter </>, and for <fr — * 1 it vanishes. Also, the heating effect 
is pronounced only for biases |A/i| > |J|, and it vanishes 
at high temperatures. For the vicinity of Ap = ± J, we 
define a characteristic temperature of a strong heating 
regime, given by Th = | J\/w (8/<p 2 _), where the function 
w(x) is defined for x > e by u>(a;) = In (x In (a;...)). Bel- 
low this temperature Th, the exponential dependence of 
exp(— J/T g) on T is replaced by a power law depen- 
dence; the temperature T competes with \Ap\ — \ J\, and 
as the latter becomes larger, the occupation probabili- 
ties cease to depend on T. For the strong heating regime 
T < Th, we find 



T 



eff 



In 1 



\J\ 



|A M |>|J|. (18) 



The current through the DD consists of an elastic and 
inelastic component, I = I c \ + T nc i. In the middle of 
the Coulomb blockade valley the current components are 
given by 



id = t- 



27 2 



(19) 



h (1 - S 2 ) 2 
e 2 7 V 

'■„„ 1 = TT— — {[©(-^ + A M) - ©(- J - A/l)] 3 Pl 

(20) 



hl-S 2 
+ [Q(J + Afi)-e(J-Afi)]p 2 } 



where 6( J) = J/(l - exp(- J/T)), and 7 = irvt L t R /E Cl 
with being the Coulomb blockade half-gap, and v 
the density of states in the leads. In the absence of heat- 
ing, when T e ff = T, the elastic component I c \ is linear 
in the applied bias Ap,, and the inelastic one 7i ne i ex- 
hibits a threshold-like switching-on at \Ap\ = \J\, for 
T < \J\ j47|. This results in steps in differential con- 
ductance versus Ap at \Ap\ — \ J\, which can be used to 
measure the singlet-triplet splitting experimentally [p5[ . 
For T <§; |J|, the step height was found to be 3 times 
larger on the singlet side than on the triplet side [|^J . In 
the strong heating regime, the differential conductance 
G = edl/dAp provides information also about the DD 
parameter 4>. We plot G(Ap) in Fig. ||, for (a) a singlet 
and (b) a triplet ground state. The dashed line shows 
the result of a calculation where the heating effects are 
neglected. We find that, in the strong heating regime, 
dlinei/dAp has a negative (positive) slope on the plateau 
\Ap\ > \J\ for the ground state being a singlet (triplet). 



The slope of the elastic component is not generic, but 
depends on the interplay between the parameters S and 
<j). However, we still find that the slope of the total con- 
ductance is not changed qualitatively by the elastic com- 
ponent over a large range of param eters calculated in the 
Hund-Mulliken method. Sec. IV A. 




°-° 45 1 2 3 

A|i/IJI 

FIG. 6. Differential conductance G = dl/dAp versus bias 
calculated for dot separation a/as = 0.7 and coupling to the 
leads given by 7 = 0.1 at a magnetic field: (a) B = and (b) 
B — 1.5 T. The solid line and the dotted line are calculated 
in the strong heating regime at T = 0.2Th and T = 0, respec- 
tively. The dashed line corresponds to neglecting the heating 
effects, i.e. T c g — > T, and was calculated for T — 0. Thus we 
see that G is monotonically increasing for J < 0, while it has 
a maximum for J > 0. 

Next, we consider the case of weakly coupled dots (4> 2 _ <C 
1) and show how the parameter <f> can be extracted from 
G(Ap) in the strong heating regime. For the singlet 
ground state, and the bias satisfying T <§; |A^i| — \ J\ <g; 
I J\l4> 2 -i we nn d 



A 



I- A 



|AH-|J| 

2\J\ 



A 



\J\ dG 



(21) 



AG dA/i' 

For the triplet ground state 
. in (|l|)^ In Fig. we plot 
2 _ in (pip at different tem- 
peratures. In the strong heating regime (T <C Th), the 



where AG = G(Ap)-G(oo) 
one should replace — ► |^ 
the rhs of the equation for 
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curve saturates at the value of <j) 2 _ with increasing bias. 
This behavior can be used to measure the interaction pa- 
rameter (j> experimentally. We note that one can reliably 
determine G(oo) experimentally as the conductance at 
Afi ^S> \J\/<p 2 _, only if the sequential tunneling processes 
can be excluded. For the middle of the Coulomb blockade 
valley, we require |J| + \A.fi\ < Ec, to avoid sequential 
tunneling via the heated excited state p6| . Furthermore, 
we note that at low temperatures Kondo resonances can 
develop at A// = ±J, invalidating Eq. (|2l]). However, 
one can avoid this problem by satisfying cither T 3> Tk 
or |A/x| — | J| ^> Tk, with Tk being the energy scale of 
the Kondo resonance. 




A|i/J 

FIG. 7. A way to measure (j) 2 _ experimentally, see Eq. 
For the calculation we used a dot separation of a fas = 1 and 
5 = (singlet ground state). Note that at T > Tk the lines 
saturate at a value which differs strongly from <fi_ . 



Kondo effect of two coupled dots at the 
singlet-triplet degeneracy point 



Attaching metallic leads to a DD can give rise to Kondo 
correlations at low temperatures. On the one hand, 
such correlations can enhance the conductance through 
a Coulomb-blockaded DD, making the transport mea- 
surements more accessible. On the other hand, Kondo 
correlations can modify the studied system and intro- 
duce a discrepancy between the measured values and the 
bare values of the system parameters. Knowledge of the 
Kondo effect in double dots would allow one to optimize 
the experimental setup, in order to obtain reliable data. 

We adopted a "poor man's" scaling approach |57]] to 
the DD system on Fig. ||and obtained an effective Hamil- 
tonian p?] ] , describing the flow (with lowering T) into the 
Kondo regime of a DD at the singlet-triplet degeneracy 
point J = 0. We found that the 4-fold degeneracy of 
the DD enhances the Kondo correlations on the Fermi 



surface, as compared to the case when the dots are de- 
tached from each other, or when the DD spins are locked 
into a spin 1 (triplet side). The Kondo temperature at 
J = is given by |Q Tk = D exp (—rj/ulo), where 
rj < 0.5 is a non-universal number dependent on the in- 
ternal features of the DD, Do ~ h~u>o is the cutoff energy 
(see Fig. |), and I = [t\ + t 2 R )/E C - We find that, at 
temperatures T < Tk, the DD undergoes a strong renor- 
malization of its energy levels, resulting in a flow of the 
exchange interaction J |47|] . Thus, we conclude that, at 
such temperatures and at J < Tk, the Kondo correla- 
tions strongly modify the coupling constant J, making 
any direct experimental measurement of the bare value 
of J hardly possible. However, at larger values of J > 0, 
one can make use of the Kondo correlations and still have 
a reliable measurement of J. 

We calculated the current through the DD at a bias 
A/i, and found the renormalization of the linear conduc- 
tance G = (dl/dA^i) | a „_q due to Kondo correlations. 
In Fig. H we plot the linear G versus the inter-dot tun- 
neling amplitude to for different values of the magnetic 
field B. At B = 0, the renormalized conductance (solid 
line) shows a sharp peak at a small value of to, owing to 
a competition between the Kondo effect of each dot with 
the adjusting lead and the antiferromagnetic exchange 
J. The peak position corresponds to J{t) ~ Tk gg| ■ At 
larger values of to, a second broader peak occurs, which 
is sensitive to applying a weak magnetic field, such that 
J > 0. We find that the broad peak is present only 
if J devia tes from At 2 H /UH by the contribution V (see 
Sec. IV A ), which comes from the long range Coulomb 
interaction (3{| (compare the solid lines with the dot- 
dashed line in Fig. ||). Note that exactly this contribu- 
tion to J is responsible for the singlet-triplet transition 
in DDs. Thus, we have shown that the long range part 
of the Coulomb interaction can be probed experimentally 
in DDs, and screening effects can be studied. 

On the left inset of Fig. ||, we plot the temperature de- 
pendence of the linear conductance calculated with tak- 
ing into account the Kondo correlations (solid line) and 
neglecting them (dotted line), the la tter c orresponding 
to the cotunneling calculation of Sec. VIA. For the case 
with Kondo correlations, we find a pronounced maximum 
in the linear G versus T at T ~ J, which can be used to 
estimate J experimentally. 

The B dependence of the linear G shows a peak at 
the singlet-triplet transition, which grows with lower- 
ing T down to T K [0, see right inset of Fig. f| Note 
that the enegry scale for the Kondo effect on the triplet 
side (J < 0) is monotonically decreasing with increas- 
ing | J | p^Jfq ] , as follows from a two stage RG procedure 
valid on that side. Furthermore, the strong coupling limit 
(not shown in Fig. ^ occurs in two stages with lowering 
T on the triplet side, resulting first in an increase and 
then, at a lower energy scale, in a decrease of the con- 
ductance pi. 
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FIG. 8. Linear conductance G at different values of B. 
Dotted lines: cotunneling contributions. Dot-dashed line: 
G vs to at B = OAT neglecting the long range part of the 
Coulomb interaction. For definiteness we keep the DD in the 
middle of the Coulomb blockade valley by adjusting the gate 
voltage V g when varying to, and choose tt = tR- Left inset: 
G vs T showing a peak at T ~ J; dotted line is the cotun- 
neling contribution. Right inset: G vs B at the singlet-triplet 
transition; the kinks in the dotted-line regions come from a 
simplified treatment of the Kondo effect crossover regions and 
will be smoothened in an exact treatment; the star denotes 
the value of B at which the singlet-triplet transition occurs 
at high temperatures (T 3> Tjf). 



VII. CONCLUSIONS 

We have described a concept for a quantum computer 
based on electron spins in quantum-confined nanostruc- 
tures, in particular quantum dots, and presented theo- 
retical proposals for manipulation, coupling and detec- 
tion of spins in such structures. We have discussed the 
requirements for coherence, switching times, read-out, 
gate operations and their actual realization. By putting 
it all together, we have illustrated how a scalable, ail- 
electronically controlled quantum computer can be envi- 
sioned. 
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